Low-frequency magnetic response of gold nanoparticles

Gold nanoparticles (AuNPs) exposed to low frequency magnetic fields have shown promise in enhancing biological processes, such as cellular reprogramming. Despite the experimental evidence, a comprehensive understanding of the underlying physical principles and the corresponding theory remains elusive. The most common hypothesis is that functionalized nanoparticles transiently amplify magnetic fields, leading to improved cellular reprogramming efficiency. However, a detailed investigation on this topic is lacking. This paper bridges this knowledge gap by conducting a comprehensive investigation on the magnetic response of surface-modified AuNPs exposed to magnetic fields with frequencies up to hundreds of MHz. Starting with the inherent properties of bulk gold material, we explore a wide range of magnetic susceptibilities that might result from the redistribution of charge carriers due to bond molecules on the particle surfaces. Through analytical models and numerical electromagnetic simulations, we examine various geometric factors that can enhance the magnetic response, including the number of particles, spatial distribution, size, and shape. Our broad investigation provides researchers with analytical and numerical estimates of the magnetic response of nanoparticles, and the associated limits that can be expected. We found that a magnetic field enhancement comparable to the incident field requires very high magnetic susceptibilities, well beyond the values measured in functionalized gold nanoparticles thus far.

Metallic nanoparticles have shown promising potential in various biomedical applications, such as drug delivery 1,2 , imaging [3][4][5] , and hyperthermia therapy [6][7][8][9] .Compared to bigger particles, nanoparticles have a high surface area to volume ratio, which makes them highly reactive to bind to other molecules 10 .Gold nanoparticles (AuNPs), specifically, have gained attention due to their unique optical 11,12 and chemical 13 properties with high biocompatibility 14,15 .In the optical frequency regime AuNPs exhibit localized surface plasmon resonances due to the collective oscillation of free electrons in response to an electromagnetic field.Furthermore, AuNPs possess antibacterial properties 16 and can catalyze chemical reactions due to their unique electronic properties 17 .For the use of electromagnetic fields (EMFs) at high frequencies (up to the visible regime) in combination with AuNPs, different physical mechanisms of interaction were already determined like radio-frequency EMFs heatinduction (as used for ablation of cancer cells) 18,19 and light-induced Raman-active vibrations (as used for the detection of cells) 20 .These and many more mechanisms of action, as heat-induction due to thermoplasmonics 21 and phothermal therapy 18,19,22 are also focuses of current research.
In the low frequency regime studies including AuNPs do usually not address the physical interaction mechanism between nanoparticles and magnetic fields [23][24][25][26][27] .However, there are hypotheses that the presence of nanoparticles causes relevant magnetic field alterations, such as changes in local field distribution or field enhancement, which lead to enhancement of biological responses [28][29][30][31] .Remarkably, application of low frequency magnetic fields with AuNPs has been shown to enhance conversion efficiency in cellular reprogramming 28 .In this regard 28 , mouse fibroblasts were cultivated on a substrate of AuNPs and treated with magnetic fields in the presence of transient reprogramming factors to mediate their direct lineage reprogramming into induced dopamine neurons.The authors found a significantly increased expression of neuronal marker genes and a 20-fold increase in conversion efficiency compared to the controls.In another study 29 AuNPs were injected in-vivo into the hippocampus of mice and stimulated high-efficient neurogenesis by exposure to magnetic fields, confirming the beneficial effect of AuNPs during magnetic field treatment.The authors in both studies 28,29 hypothesize that significant transient magnetization of the particles is causing the biological reactions.This theory is also taken up particle, which will serve as a foundation for our subsequent analyses.Also as a basis for subsequent investigations, we show the magnetic response of an individual AuNP with bulk gold material properties and investigate the effect of the particle size on the magnetic response.Finally, we consider the interaction of two spherical particles as a function of their magnetic susceptibilities ( −0.9 ≤ χ p ≤ 10 ) as this provides the basis of both: the investigation of interaction of different shapes of particles and the investigation of interaction of several particles.

Particle polarizabilities
The magnetic dipole moment vector m of a spherical particle which is exposed to an incident magnetic field H inc can be expressed with two magnetic polarizability contributions: α m and α e .α m is the magnetic polarizability due to magnetization of the particle, and α e is the magnetic polarizability due to conductive and displacement currents within the particle (details on the calculations can be found in the methods section): The experimental studies described in the introduction were conducted at frequencies f < 1MHz [23][24][25][26][27][28][29][30][31] .Since we will consider nanoparticles with sizes not exceeding several tens of nanometers, the wavelength of the external field is much bigger than the considered particles and R p k 0 ≪ 1 , where R p is the particle radius and k 0 is the wavenumber in free space.Under these conditions, for calculation of α m , the quasi-static approximation 44 can be applied leading to the following expression 45 where χ p is the volume magnetic susceptibility with the relative magnetic permeability being µ r = 1 + χ p and V = 4πR 3 p /3 is the volume of the particle.In the quasi-static approximation (negligible retardation), with R p k 0 ≪ 1 , α e can be expressed with 46,47 k = k 0 √ ε p is the wavenumber in the magnetized volume V, and ε p = ε p,r + iσ p /(2πf ε 0 ) is the relative permittivity of the particle accounting the dielectric part contribution ε p,r and the specific conductivity σ p .ε 0 is the vacuum dielectric constant, i is the imaginary unit.Equations ( 2) and ( 3), show that while α m is independ- ent of frequency, α e is frequency-dependent due to the wave number dependence.The ratio between the two polarizabilities can be presented as Since α e /α m ∼ R 2 p k 2 0 , it can be expected that in the low-frequency range the main mechanism of interaction of gold particles with an external electromagnetic field is their quasi-static magnetization.Figure 1 illustrates the ratio of the magnitudes of the two magnetic polarizability contributions, |α e |/|α m | , as a function of frequency f and particle radius R p .To model the behavior of a spherical gold particle, we used the material parameters for bulk gold, specifically the electrical conductivity σ p = 41 MS/m , the relative permittivity ε r,p = 1 , and the magnetic susceptibility χ p = −3.4× 10 −5 .Figure 1 shows that for R p ≤ 100 nm and f ≤ 100 MHz , |α m | ≥ |α e | .
Figure 1.Size-and frequency dependent ratio of the two magnetic polarizability contributions in an AuNP with bulk gold material properties.
Since the experimental studies described in the introduction were conducted at frequencies f < 1MHz [23][24][25][26][27][28][29][30][31] , we consider in this article that |α m | ≫ |α e | , and the magnetic dipole moment induced in the nanoparticle can be obtained with the static approximation m = α m H inc , neglecting a contribution of the dynamic response.Note that the analytically obtained result presented in Fig. 1 has been validated with numerical electromagnetic simulations in Comsol Multiphysics (details on the simulations are presented in the "Methods" section).

Response of a single particle with bulk gold material properties
In the static approximation, the expression for the magnetic field outside a spherical particle with magnetic dipole moment m is where r represents the vector connecting the center of the sphere and the point at which the magnetic field is being evaluated.With (5), the following equations can be obtained When considering this formulation from (5) it becomes clear that depending on the magnetic susceptibility χ p and the angle of r relative to H inc , the magnetic field H can be enhanced or lowered compared to the incident field [note the different signs in the Eq. ( 5)].
The magnetic field inside the sphere can be calculated with 45 Results in Fig. 2 were calculated for an AuNP with radius R p = 10 nm that is exposed to H inc directed along the z-axis of the Cartesian coordinate system with the origin at the particle center.Figure 2a illustrates the normalized change of the magnitude of the total magnetic field with respect to the incident field: being the total magnetic field magnitude).By applying this definition, the magnitude of H norm shows the relative change in H compared to H inc due to the particle and the sign of H norm indicates whether H is lower or higher than H inc .Within this article H norm is referred to as the secondary mag- netic field.In Fig. 2a the diamagnetic AuNP shows maximum enhancement of the secondary magnetic field where the incident magnetic field is oriented tangential to the surface of the sphere, and the secondary magnetic field is weakest where the incident magnetic field is oriented perpendicularly relative to the surface of the sphere.The inverse would be the case for a particle with positive magnetic susceptibility.However, the maximum magnetic field enhancement is in the order of 10 −5 relative to the incident field.Figure 2b shows the gradient of Fig. 2a outside the particle, where the gradient magnetic field reaches the maximum order of 10 4 .This relatively high magnetic field gradient is confined to the immediate vicinity of the sphere.When considering (2) the polarizability of the particle depends on its radius R p and its susceptibility χ p . Figure 2c shows the secondary magnetic field distribution of the same arrangement as in Fig. 2a along the x-axis for different particle sizes.It can be seen that the maximum magnetic field remains the same for all three particle sizes while the maximum gradient of the magnetic field varies significantly with the size of the particle.The maximum magnetic field gradient decreases with increasing particle size.
From (2) it can be deduced that for the spherical particle the magnitude of the polarizability becomes bigger with increasing magnitude |χ p | .While for |χ p | ≪ 3 the polarizability can be approximated with α m = V χ p , saturation is reached for |χ p | ≫ 3 ( α m = 3V ).While the quantitative magnetic field distribution changes for different magnetic susceptibilities, the qualitative magnetic field distribution around a spherical particle (location of minima and maxima, and overall field distribution) remains similar to the results shown.

Coupling of two spherical particles
Since in biological applications the surface of a nanoparticle is usually modified to improve its biocompatibility, targetability, stability, and functionality 48 , the magnetic susceptibility of the particle can differ significantly from the susceptibility of bulk gold material due to charge carrier redistribution.Diamagnets obtain a magnetic susceptibility smaller than zero, while superconductors are considered ideal diamagnets with a magnetic susceptibility of − 1 49,50 .Paramagnetic materials obtain magnetic susceptibilities higher than zero and the magnetic susceptibility of ferromagnetic materials can be in the range of hundreds of thousands as in the case of iron 51 .We choose a wide range of magnetic susceptibilities from −0.9 ≤ χ p ≤ 10 for our investigations, as to the best of our knowledge, the possible limits of magnetic susceptibilities for surface modified AuNPs are not yet determined 41 .This range is probably wider than realistic magnetic susceptibility values for modified AuNPs; however, it allows us to investigate the limits of possible magnetic field interaction.
In the following, the interaction of two spherical particles as a function of their magnetic susceptibility is analysed.
Figure 3a-d show two touching AuNPs arranged along the y-axis.Furthermore, in Fig. 3a and b, the incident magnetic field is oriented in y-direction.Applying the general Eqs.(30) and (31) from the methods section, the magnetic moment of each particle is oriented in y-direction and its magnitude can be expressed as while ( 8) H 2,y is the additional magnetic field contribution from the second particle at the position of the first particle and m 2,y is the magnitude of magnetic dipole moment of the second particle.Since the particles are identical, m 1,y = m 2,y .Thus, by considering the polarizability from (2) and reformulation of (8) the following equation can be obtained Thus, the magnitude of the magnetic dipole moment of each particle from Fig. 3a and b can be calculated with the superposition of the following contributions The index zero refers to the dipole moment of each particle in the absence of interaction with the other particle.The index ' coupling' indicates the additional magnetic dipole moment resulting from the interaction between the particles.
Considering Fig. 3c and d, the particles are arranged along the y-axis and the incident magnetic field is oriented in z-direction.By application of the same approach as before, the magnitude of the magnetic dipole moment of each particle from Fig. 3c and d

can be calculated with the superposition of the following contributions
The applicability of the presented formulas was validated by comparison with numerically obtained magnetic dipole moments.The comparison is shown in the Methods section, where it can be seen that the formulas result in very good agreement with numerically obtained magnetic dipole moments.
In Fig. 3a and b the ratio m y,coupling /m y,0 is plotted and in Figure 3c and d the ratio m z,coupling /m z,0 is shown.The respective ratios can be positive or negative, depending on the susceptibility χ p and the orientation of the magnetic field vector relative to the particle arrangement.This means that the interaction of the particles can cause an increase or a decrease of magnetic dipole moment.Additionally, in all cases, the coupling increases with increasing magnitude |χ p | , while in the case of para-and ferromagnetic particles the onset of saturation is visible around χ p = 10.
It should be noted that for higher particle distances the coupling of the particles becomes weaker until the magnetic dipole moment can be approximated with the magnetic dipole moment of individual particles ( m = α m H inc ).This is further shown in the "Methods" section.

Prolate and oblate surface-modified particles
In the following, the interaction of elliptical particles will be treated and compared to the spherical case.Ellipsoidal AuNPs are among the particle shapes commonly used in the frame of biomedicine 20 .Furthermore, it was found that ellipsoidal particles thiolated with polyethylene glycol can exhibit strong diamagnetic properties 37 .We fist concentrate on the magnetic response of individual prolate elliptical particles.We investigate the response as a function of the different axis ratios as well as a function of different magnetic susceptibilities.Subsequently, we show the results for oblate elliptical particles.The investigations of the individual particles allow a target-oriented investigation of the coupling of two particles which are spherical, prolate or oblate elliptical.Since we like to explore the limits of magnetic response, the target is to achieve a maximum in magnetization of the particles.Therefore, we investigate the magnetization as a function of two interacting particles with the different shapes.

Prolate elliptical particles
Figure 4 depicts prolate elliptical particles with semi-axes of lengths l x , l y and l z , where l z > l y and l y = l x .The indices denote the orientation of the length parallel to the respective axis, while l z /l y = p .The parameter p equals integer values in the range 2 ≤ p ≤ 100.
Figure 4 shows the prolate elliptical particle with the incident magnetic field vector oriented either along the major axis of the particle (Fig. 4c,d) or one of the minor axes (Fig. 4a,b).The magnitude of magnetization M i of one particle can be calculated with 52 where N i is a shape dependent (demagnetization) factor, since the shape of the particle affects the magnetic field distribution.The Cartesian axes are elements of the index i and indicate the orientation of the magnetization.Since a sphere has point-symmetry relative to its center, its demagnetization factor is constant in all directions (9) and equals N i = 1/3 53 .The formulas ( 35) and (36) in the "Methods" section allow the calculation of the demag- netization factors for prolate elliptical particles.The graphs in Fig. 4 show the relationship between the respective magnetization of the particle relative to the magnetization of a sphere, as a function of the particle's susceptibility, and different ratios of the major axis length to the minor axis length ( l z /l y = p ). Depending on the magnetic susceptibility and orientation of the magnetic field vector relative to the particle, the magnetization of the prolate particle can be greater or lower than the magnetization of a sphere with the same susceptibility.The biggest enhancement in magnetization is found for χ p > 1 with the magnetic field vector oriented parallel to the particle's major axis (Fig. 4d).

Oblate elliptical particles
Similarly to the results from Fig. 4, the graphs in Fig. 5 show the magnetization of oblate particles relative to the magnetization of a sphere with the same magnetic susceptibility.The magnetization of the oblate elliptical particle is calculated with (13) and by application of the formulas (37) and (38) for the demagnetization factor.The highest magnetization occurs when the incident magnetic field vector is perpendicular to the minor axis of the oblate particle and χ p > 1 (Fig. 5d).
When comparing Fig. 5d to 4d, it is evident that the prolate particle induces a much higher magnetization at χ p > 1 than the oblate particle.Interestingly, when comparing the diamagnetic particles in Fig. 5a and 4a, the oblate particle yields a higher magnetization.

Coupling of the differently shaped particles
The comparison of Figs. 4 and 5 allows assumptions of the magnetic field interaction of the individual particles, especially the quantification of the impact of the three different shapes (spherical, prolate and oblate ellipsoidal).However, the interaction of elliptical particles is yet to be investigated and compared to the interaction of spherical particles.
In this subsection, we investigate the limits of magnetic responses and search for a configuration of two identical particles that results in maximum magnetization.First, we consider diamagnetic particles.Figure 3c shows the configuration that generates the highest magnetic dipole moment for diamagnetic spherical particles.Positive coupling of diamagnetic particles can be achieved when the incident magnetic field is perpendicular to the line connecting the particle centers (e.g. the incident magnetic field in z-direction and the arrangement of particles in y-direction as in Fig. 3c).Additionally, for elliptical diamagnetic particles, magnetization is highest when the magnetic field is oriented along one of the semi axes (as seen in Figs.4a and 5a).The maximum coupling of particles is achieved when they are brought as closely together as possible.Based on these information, Fig. 6 shows the configurations of diamagnetic particles with different shapes that result in maximum magnetization.The graph illustrates M y,norm , which is the magnitude of magnetization in y-direction, normalized to the incident magnetic field and can be expressed by application of the magnetic dipole moment: Note that in this definition the magnetization is normalized to the incident magnetic field to avoid the dependence of the parameter from the incident field.When approximating the particles as point dipoles, the general Eqs. ( 30) and (31) from the methods section can be applied to calculate the dipole moments.This results in a similar approach as in the case of the coupling of two spheres and the magnetic dipole moments of the particles from Fig. 6 can be expressed with α m,y = χ p V /(1 + N y χ p ) is the polarizability of the particle in y-direction.Solving (15) for m y and inserting the result into ( 14), the following equation can be obtained In Fig. 6, the different colors represent the different susceptibilities.M y,norm is plotted as a function of the configurational parameter n, which represents the ratio of the lengths of the axes The numerically obtained results from Comsol Multiphysics are depicted using dashed lines, while the results obtained analytically with ( 16) are represented with solid lines.Generally, the numerically and analytically obtained results are in good agreement, with the magnetization increasing with higher magnetic susceptibility magnitudes, which is unsurprising.It can be noticed, that for χ p ≤ −10 −1 the analytical results are obtained for n > 0.3 .The reason thereof is that a limit for the applicability of ( 31) is reached, which is used for the derivation of (16).For a very long prolate particle's short-range area, the point dipole approximation is not valid.However, for all other susceptibility values and shapes, there is significant agreement between the analytical and numerical results.
At χ p = −0.9, it is noticeable that the oblate particle with the maximum ratio of major axis to minor axis (n) generates the highest magnetization.This trend can also be observed in magnified inset on the other curves.Additionally, in magnified inset, it is observed that both prolate and oblate particles yield higher magnetization than spherical particles.
Figure 7 illustrates a result similar to that shown in Fig. 6, but for the coupling of para-/ferromagnetic particles.The principle mathematical approach is similar to the derivation of ( 16) and the normalized magnetization is Figure 7 demonstrates excellent agreement between the numerically and analytically obtained results.The prolate particles exhibit the highest magnetization which is in contrast to the diamagnetic particles, where the oblate particles yield to highest magnetization.

Example: spherical arrays of particles
As already mentioned in the introduction one major study used surface modified AuNPs for biological reprogramming 28 .There, fibroblast cells were planted on an array of spherical particles and exposed to a magnetic field 28 .The arrangement of particles raises the question: how do the particles interact with each other?Specifically, we seek to determine the number of particles that interact with one another in such an array.To determine the number of particles that interact with each other in an array, we first focus on a line of particles, as illustrated in Fig. 8.To enhance the magnetization of the diamagnetic particles in the array due to their interactions, as in the result of Fig. 3c, we applied a magnetic field perpendicular to the direct connection vectors between the particle centers.We gradually increased the number of particles in a line from N = 2 to N = 6 , as shown in Fig. 8.The figure shows the ratio of the maximum magnetization ( M y,array ) within the line of particles to the magnetization of a single particle without any interactions ( M y,0 ).The asterisks indicate the corresponding M y,array /M y,0 values for each particle number, while the dashed lines show how this ratio changes with N. The colors in the figure represent different magnetic susceptibility values.The susceptibility curves for χ p ≥ −10 −2 (16) x + α m,y .
( www.nature.com/scientificreports/overlap, indicating a lack of significant increase in magnetization due to particle interaction at these levels.However, for susceptibility values lower than χ p ≤ −10 −1 , there is a relevant increase in M y,array /M y,0 , showcasing meaningful particle interaction.The maximum interaction between particles saturates at N = 3 for χ p = −10 −1 , whilst for χ p = −0.9, saturation is reached at N = 5 particles.Based on the susceptibility values considered, it can be inferred that a maximum of five particles aligned in a single row can result in relevant enhancement in magnetization.Expanding the number of particles in the row only leads to an increase in the area covered by the maximum magnetization, but does not enhance the maximum magnetization itself.Similar to Figs. 8 and 9 shows results for particle interaction in the form of M y,array /M y,0 for para-/ ferro- magnetic particles.The arrangement of particles relative to the incident magnetic field is chosen on the basis of the results from Fig. 3a to obtain an enhancement of magnetization.The results show that the onset of increased M y,array /M y,0 begins at χ p ≥ 10 −1 , and the number of interacting particles increases with higher magnetic sus- ceptibility values.Saturation is reached at N = 9 for χ p = 10.
Returning to an array of diamagnetic particles as used in a study for biological reprogramming of cells 28 , we will now investigate the limits of magnetic field alteration due to the particles interaction.Figure 10a shows an array of AuNPs, six particles wide and six particles tall, exposed to a magnetic field perpendicular to the array plane.The particles have a diameter of R p = 10 nm and are placed directly beneath each other without any addi- tional spacing.The magnetic field H norm is determined within the black rectangular area displayed in Fig. 10a, which is then illustrated in Fig. 10b.The same arrangement is analyzed for an array of 10 by 10 nanoparticles, as depicted in Fig. 10c.Note that Fig. 10b and c use the same scale.Comparing the two plots, it is evident that there is no significant enhancement in the magnetic field due to more particles.Although the alignment of the particles enables an enhancement in magnetization, it leads to an increase in magnetic field between the particles and a reduction in magnetic field on the particles' surface.This results in a range of −0.5 ≤ H norm ≤ 0 , showing that the magnetic field above the array is lower than the incoming field.
The magnetic field gradient of the arrangement from Fig. 10a and b is shown in Fig. 10d.The magnetic field gradient peaks at approximately 10 8 m −1 .
The magnetic field distribution from Fig. 2a shows that the secondary magnetic field around an AuNP is maximum at the positions where the incident magnetic field vector is only tangential to the surface of the particle.To obtain a magnetic field enhancement above an array of diamagnetic particles, this characteristic is exploited.In contrast to the arrangement from Fig. 10a and 11a shows a configuration where a partial enhancement of the magnetic field on top of the particles is excited.The incoming magnetic field is running parallel to the z-direction.The particles are closely positioned in the x-direction while there is a 60 nm distance between neighboring particles in the z-direction.The spacing results in a decrease of negative coupling (as in the result of Fig. 3a) between the diamagnetic particles.By implementing this configuration, the magnetic field is increased by a factor of 0.3 directly above the particles.
The findings from Fig. 6 indicate that the magnetization and magnetic field interaction of particles could be enhanced by utilizing prolate/oblate particles with low/ high configurational parameter n.However, when dealing with spherical diamagnetic particles, achieving magnetic field enhancement in the range of the incident magnetic field necessitates magnetic susceptibility values in the range of χ p = −0.9 and a very specific particle arrangement, as chosen in Fig. 11.

Discussion
We presented analytical formulas for the calculation of the magnetic response of individual spherical particles for frequencies up to several hundreds of Megahertz.These formulas provided the basis of our following investigations.All our analytically obtained results are validated with numerical electromagnetic simulations in Comsol Multiphysics.We showed that the magnetization of an AuNP is dominated by the polarizability α m as long as the magnetic field excitation fulfills f < 100 MHz and R p ≤ 100 nm .We showed the magnetic field distribution within and around an AuNP with bulk gold material properties, as well as its magnetic field gradient.These results showed that the magnetic field around an AuNP is partially enhanced and partially decreased, depending on the considered position.However, the alteration in the magnetic field is approximately five orders of magnitude smaller than the incident magnetic field while the maximum magnetic field gradient is in the order of 10 5 .Considering that the magnetic energy E of a particle is 54 (µ 0 is the vacuum permeability) the magnetic field contribution of an AuNP with bulk gold material parameters is too small to cause significant energy alterations.Furthermore, the magnetic force per particle volume f is 54,55 Thus, to achieve a magnetic force density in the range of gravitational force density which is approximately f g = 10 4 N/m 3 55 , very high magnetic fields in combination with magnetic field gradients are necessary.However, we also showed that the gradient of the magnetic field can be tuned with the particle size.
Furthermore, the magnetic dipole moments of two neighboring spherical particles in contact can be enhanced or lowered due to their interaction.In this regard, the determining factors are orientation of the magnetic field relative to the arrangement of nanoparticles, and the susceptibility of the particles.
Since the magnetization of a particle is proportional to its magnetic field response, we investigated the magnetization of prolate and oblate ellipsoids relative to the magnetization of a spherical particle.We considered different axis ratios for the elliptical shapes and also accounted for different magnetic susceptibilities.Thereby, we showed that the maximum magnetization (hence maximum magnetic field) can be reached for prolate para-/ ferromagnetic ellipsoids and oblate diamagnetic ellipsoids.Furthermore, depending on the orientation of the magnetic field relative to the elliptical particle, the magnetization can also be lower than the magnetization of spheres.
Subsequently, we investigated and compared the magnetic coupling of two particles with the different shapes.Here, we saw that in the cases of positive interaction (interaction causing an increase in magnetization), the elliptical particles show stronger interaction than spherical particles.We showed that for diamagnetic particles, the oblate ellipsoids show the strongest interaction while for para-/ferromagnetic particles the prolate ellipsoids interact most strongly.
Emanating from an example taken from literature, we investigated the interaction of spherical particles in an array.We found that the determining coupling factor is the magnitude of the magnetic susceptibility: the higher the magnitude |χ p | , the more particles can interact with each other.However, for the considered susceptibility range −0.9 ≤ χ p ≤ 10 in case of diamagnetic particles maximum five particles can couple effectively in one row and in case of ferromagnetic particles, maximum fifteen particles can couple effectively in one row.Finally, we exemplary investigated the magnetic field and the magnetic field gradient above arrays of diamagnetic nanoparticles.The main findings of these investigations are that for the case of positive coupling of particles, depending on the orientation of the incident magnetic field relative to the array surface, the incident magnetic field can be enhanced or lowered.Our results show that an extreme diamagnetic susceptibility of χ p = −0.9can achieve a magnetic field enhancement in the range of the incident magnetic field.Since there is the hypothesis that a relevant increase in magnetic field was reached with strong diamagnetic AuNPs experimentally 28 , we determine that the magnetic susceptibility for such an increase in the magnetic field must be in the range of χ p = −0.9 .As mentioned before, the limit of possible susceptibility values for surface-modified AuNPs is yet to be determined.However, a value in the range of χ p = −0.9 is beyond realistic susceptibility values found in literature thus far (see Introduction).
We would like to emphasize that to the best of our knowledge, the exact physics of the magnetism in functionalized AuNPs is not finally understood.Our study provides researchers with a methodology to estimate the magnetic response of functionalized AuNPs with known magnetic susceptibility.We identified formulas which are applicable, and showed their validity with numerical electromagnetic simulations.In contrast to numerical methods, analytical approaches directly show the connection between different parameters (as shape, size and number of particles) and the magnetic response.The understanding of the low-frequency magnetic response of AuNPs is essential for further determination of the interaction of the magnetic field energy with biological cells and cell compartments.Our findings contribute to a better understanding of the limits, characteristics, and impact of AuNPs on the magnetic field response, with potential applications in biomedical research and treatment.Therefore, our investigation provides important guidance and contributes to a better understanding of the low-frequency magnetic response of AuNPs.
However, we believe that further research is needed to comprehend the physics behind the enhanced reprogramming efficiency observed 28 when using functionalized AuNPs in combination with magnetic fields.The original hypothesis that an enhanced magnetic field due to AuNPs causes a higher reprogramming efficiency seems unlikely, because for realistic magnetic susceptibility values the enhancement of the magnetic field is orders below the incident magnetic field.Since the exact physics behind functionalized AuNPs are not known yet, a further direction to be investigated is the fundamental investigations on the properties of materials, especially the exact relation between factors that affect the magnetic susceptibility and the magnetic susceptibility itself.Furthermore, the relevance of non-uniform magnetic field distributions for the interaction of AuNPs and cells needs to be investigated since AuNPs cause relatively high magnetic field gradients in their close vicinity.Finally, (19

Numerical electromagnetic simulations
Numerical simulations were performed using the Magnetic Fields (mf) interface under the AC/DC Electromagnetic Fields branch in Comsol Multiphysics (COMSOL AB, Sweden, https:// www.comsol.com).Due to high water content, the effective magnetic susceptibility of cells, including major compartments like the cytosol, is in the range of 10 −655,56 .Lipid structures in cells have even lower magnetic susceptibilities than compartments with high water content 55 .However, due to their higher magnetic susceptibility, the magnetic response of AuNPs to low-frequency magnetic fields is generally much stronger than the response of eukaryotic cells.Therefore, in our investigations, we neglect the specific cellular environment and consider AuNPs located in air.The geometry consisted of the gold nanoparticles with the corresponding shapes and material properties for the respective investigation, located inside an (spherical) environment with vacuum material properties.An additional layer around the environment built up the infinite element domain to absorb outgoing EMFs.The incident magnetic field was excited by setting a background magnetic vector potential which resulted in the target incident magnetic field.However, when evaluating the results, the total magnetic field was normalized to the incident magnetic field.The target frequencies of the respective investigations were set in the study section and material sweeps were conducted where necessary.Magnetic dipole moments were determined numerically in the derived values section, by volume integration of the magnetization.

Fundamental analytical calculations
Ampère's circuital law and the magnetic polarization definition Ampère's circuital law in differential form and in frequency domain describes the rotation of the magnetic field H as the sum of different current density contributions, In Eq. ( 21) E is the electric field vector, σ is the electrical conductivity and ε = ε 0 ε p,r is the material permit- tivity, while M is the magnetization of it.The first summand in (21) describes the conductive and displacement current density j e and the second summand describes a current density j m due to magnetization.The magnetic dipole moment m of the considered volume can be expressed with a volume integral over the current densities from (21), V is the volume of the particle, and r is the position vector.Furthermore, the magnetic dipole moment is generally defined as the product of a magnetic polarizabilitiy and the incident magnetic field.We define α m as the magnetic polarizability which is connected to the magnetization M and α e as the magnetic polarizability which is connected to the conductive and displacement currents.
In order to obtain expressions for polarizabilities, it is convenient to consider separately problems with different current sources.In the case of j m = ∇ × M , we can make the following transformations 44,45 where, V = 4πR 3 /3 (R is the particle radius), and we used the solution to the static problem of a magnetic sphere with relative magnetic permeability µ r in a constant external magnetic field H inc 45 .Note that this solution can be used because we assume from the very beginning that the particle is much smaller than the wavelength of the external field.Thus, we can write in this approximation that In general the magnetic dipole moment contribution due to magnetization M has a static contribution and a frequency-dependent contribution 44 .However, in the considered frequency regime from Fig. 1, the static contribution is dominant leading to the static magnetic polarizability α m from (2).
The case of j e = (σ − iωε)E is more complicate and can not be considered in the static approximation because we must take into account the connection between electric and magnetic fields corresponding to the Maxwell equation: where k = ω √ ε 0 µ 0 ε p , and ε p = ε p,r + iσ/(ωε 0 ) .Finally, for the magnetic field inside and outside the particle we have the equations where to write the equation for the inside we used that ∇ × ∇× = −� + ∇(∇•) and ∇ • H = 0 .The solution of (27) for a spherical particle has been considered elsewhere 47 , where it is shown that the total magnetic field outside a particle can be presented as a superposition of the incident external field and the field generated by the particle's magnetic dipole moment m e determined by the polarizability If the particle size is so small that |Rk| ≪ 1 then This coincides with (3).

Main formulas for interacting particles in the magnetic dipole coupling model
The magnetic dipole moment of one individual particle is excited by the incident magnetic field.When other particles are added, they change the initial magnetic field which is experienced by the first particle and influence its magnetic dipole moment.In a configuration of N identical particles, where N ∈ N, N ≥ 2 , the magnetic moment of each particle can be calculated with The index n indicates the n-th particle and the index k indicates the other particles which interact with it.The magnetic field contribution of the k-th particle can be calculated with 57 where r = |r|.

Numerical validation of the formulas for the magnetic dipole moment of two touching spherical particles
Figure 12a shows two touching diamagnetic particles which are arranged along the y-axis and which are exposed to an incident magnetic field parallel to the z-axis.This configuration equals the arrangement of the particles from Fig. 3c, which resulted in positive interaction of particles.Due to symmetry, the magnetic dipole moment vector m of both particles is equal in magnitude and is directed in negative z-direction.To validate the ana- lytical formulas derived for the determination of the magnetic dipole moment, the analytically obtained magnetic dipole moment |m anal | was compared to the numerically obtained magnitude |m num | .In Fig. 12a  Figure 12b shows a similar graph but for two para-/ ferromagnetic particles in the same arrangement as shown in Fig. 3b. Figure 12b shows that for susceptibilities up to χ p ≈ 1 , the analytically and numerically obtained results ∇ × H = 0, and ∇ • H = 0 (outside)  The two graphs from Fig. 12 show that in general there is good agreement between the analytically obtained magnetic dipole moment and the numerically obtained magnetic moment, which speaks in favor for the presented analytical formulas.
Furthermore, Fig. 13 shows the numerical EM simulation of diamagnetic spherical particles with χ p = −0.9 at different distances d.Additionally, Fig. 13 shows the numerical EM simulation results of point dipoles approximating the spherical particles at respective distances with a magnetic dipole moment m .All plots from Fig. 13 show H norm in the same scale.The magnetic field distributions resulting from the point dipoles approximation and the respective magnetic field distributions resulting from the spherical particles are in very good agreement when considering the regions outside the sphere.In the arrangements of particles from Fig. 13 the maximum magnetic field enhancement can be obtained when the particles are brought as closely together as possible ( d = 20 nm ).When the distance between the particles is increased, the magnetic field enhancement becomes smaller until it equals the maximum magnetic field enhancement of individual particles.This can also be shown analytically.Considering (30) with (31) the magnetization m of one point dipole from Fig. 13 can be expressed with With Eq. ( 32), the magnitude m can be obtained as By rearrangement of ( 33), the ratio of the magnetic dipole moment m and the magnetic dipole moment of a single particle in absence of the other particle m 0 = |α m H inc | can be expressed as a function of the distance d. Figure 14a shows the ratio m/m 0 obtained analytically with (33) as a function of the distance d.   the numerically obtained solution of the ratio is plotted for the three distances from Fig. 13. Figure 14b shows H norm at the origin of the coordinate system from Fig. 13 obtained analytically with Eq. ( 5) and with numerical simulations at the distances from Fig. 13.It can be seen that the magnitude of H norm correlates with the ratio m/m 0 .Furthermore, there is good agreement between the numerical and analytical results.

Demagnetization formulas for prolate elliptical particles
The demagnetization factors of a prolate elliptical particles in the three Cartesian directions can be calculated with 58

Demagnetization formulas for oblate elliptical particles
The demagnetization factors of the oblate elliptical particles can be calculated with 58

Figure 2 .
Figure 2. (a) Secondary magnetic field H norm = (|H| − |H inc |)/|H inc | and (b) its gradient.Both plots in the xz-plane.(c) H norm along the x-axis for different AuNP sizes.All plots for bulk gold material properties.

Figure 3 .
Figure 3. Magnetic dipole moment contribution due to coupling of two identical spherical AuNPs normalized to the magnetic dipole moment of a single particle.(a) Diamagnetic particles with incident magnetic field in y-direction.(b) Para-/ferromagnetic particles with incident magnetic field in y-direction.(c) Diamagnetic particles with incident magnetic field in z-direction.(d) Para-/ferromagnetic particles with incident magnetic field in z-direction.

( 14 )Figure 4 .
Figure 4. Relative magnetization in a prolate AuNP.(a) Diamagnetic particle with incident magnetic field in y-direction.(b) Para-/ferromagnetic particle with incident magnetic field in y-direction.(c) Diamagnetic particle with incident magnetic field in z-direction.(d) Para-/ferromagnetic particle with incident magnetic field in z-direction.

Figure 5 .Figure 6 .
Figure 5. Relative magnetization in an oblate AuNP.(a) Diamagnetic particle with incident magnetic field in y-direction.(b) Para-/ferromagnetic particle with incident magnetic field in y-direction.(c) Diamagnetic particle with incident magnetic field in z-direction.(d) Para-/ferromagnetic particle with incident magnetic field in z-direction.

Figure 7 .
Figure 7. Normalized magnetization of para-/ ferromagnetic particles with different shapes and axis ratios.Solid lines are results obtained with presented analytical formulas.Dashed lines are results obtained with numerical EM simulations.

Figure 8 .Figure 9 .
Figure 8. Maximum magnetization in a line of spherical diamagnetic particles normalized to the magnetization of a single particle without in absence of the other particles.

Figure 10 .
Figure 10.Array of spherical diamagnetic particles.(a) Schematic illustration of the 6 × 6 particles array.(b) Secondary magnetic field distribution above the array.(c) Secondary magnetic field distribution above a 10 × 10 particles array.(d) Gradient magnetic field of the result from (b).

Figure 11 .
Figure 11.Modified array of spherical diamagnetic particles.(a) Schematic illustration.(b) Secondary magnetic field distribution above the array.

Figure 12 .
Figure 12.Ratio of analytically obtained magnetic dipole magnitude |m anal | to numerically obtained magnetic dipole moment magnitude |m num | as a function of magnetic susceptibility. (a) Interaction of particles according to Fig. 3c and (b) interaction of particles according to Fig. 3b.

Figure 13 .
Figure 13.Numerical EM simulation of spherical particles at different distances d and comparison with numerical simulation of magnetic point dipoles with the same dipole moment m as the spherical particles at the same distance.View on the longitudinal axis cross-section of the sphere through the center ( z = 0).

Figure 14 .
Figure 14.Numerical EM simulation of spherical particles at different distances d and comparison with numerical simulation of magnetic point dipoles with the same dipole moment m as the spherical particles at the same distance.View on the longitudinal axis cross-section of the sphere through the center ( z = 0).
f = µ 0 χ p (H • ∇)H.to understand the interaction of magnetic fields with surface-modified AuNPs and cells, different modelling principles considering multiphysics can be applied, e.g.molecular dynamics or density functional theory.